The Heisenberg antiferromagnet on the kagome lattice with arbitrary spin: 
A high-order coupled cluster treatment 
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Starting with the \/3 x y^3 and the q = states as reference states we use the coupled cluster 
method to high orders of approximation to investigate the ground state of the Heisenberg antiferro- 
magnet on the kagome lattice for spin quantum numbers s = 1/2, 1,3/2, 2, 5/2, and 3. Our data for 
the ground-state energy for s = 1/2 are in good agreement with recent large-scale density-matrix 
renormalization group and exact diagonalization data. We find that the ground-state selection de- 
pends on the spin quantum number s. While for the extreme quantum case, s — 1/2, the q — 
state is energetically favored by quantum fluctuations, for any s > 1/2 the \/3x \/Z state is selected. 
For both the %/3 x v3 and the q — states the magnetic order is strongly suppressed by quantum 
fluctuations. Within our coupled cluster method we get vanishing values for the order parameter 
(sublattice magnetization) M for s = 1/2 and s = 1, but (small) nonzero values for M for s > 1. 
Using the data for the ground-state energy and the order parameter for s = 3/2,2,5/2, and 3 we 
also estimate the leading quantum corrections to the classical values. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.50.Ee 



I. INTRODUCTION 



The investigation of the low-energy physics of the 
Heisenberg antiferromagnet (HAFM) 



H = 



E 1 

<i,j> 



(1) 



on the kagome lattice is one of the most challenging prob- 
lems in the field of frustrated quantum magnetism. The 
sum over runs over all nearest-neighbor pairs of 

sites on the lattice, counting each bond once only, and 
Si = (sf, s?, sf ) is the spin operator on site i. Although, 
there has been an intensive discussion of the problem 
over many years appl ying various theoretical methods, 
(and see, e.g., Refs. lfl-1341). no conclusive answer on the 
nature of the ground state (GS) and the existence of a 
spin gap has been found. 

While for many years a spin-liquid GS was 
favored ) 5 ' 10 ' 12 recently arguments have been given for 
a valence-bond crystal GS with a large unit cell of 36 
sites that breaks the symmetry of the underlying kagome 
latticej 16 ' 20 ' 21 However, very recently this valence-bond 
picture has been rechecked by large-scale numeric a 28 ' 29 ' 31 
and once again the spin-liquid GS is favored. 

Although, large-scale density-matrix renormalization 
group (DMRG) and exact diagonalization (ED) calcula- 
tions seem to be most effective to study the low-energy 
physics of the kagome HAFM, complementary methods 
(and see, e.g., Refs.HglalHal, are highly desirable to 
shed further light on the challenging problem. 

A method which has been successfully applied to 
strongly frustrated quantum magnets is the coupled clus- 
ter method (CCM) (and see, e.g., Refs. H^il)!©. In 
the present paper we apply the CCM in high orders of 
approximation to the kagome HAFM. 




FIG. 1: Illustration of the ^3 x ^3 (left) and the q = 
(right) classical GS of the kagome HAFM. 



II. COUPLED CLUSTER METHOD 

For the sake of brevity we illustrate here only some 
relevant features of the coupled cluster method (CCM). 
For more general information on the methodology of the 
CCM, see, e.g., Refs. Sill 

We first mention that the CCM approach yields results 
directly in the thermodynamic limit N — > oo, where N 
is the number of lattice sites (and hence spins). The 
starting point for a CCM calculation is the choice of a 
normalized reference state |$) that is typically a classical 
GS of the model. For the kagome HAFM we choose the 
\/3 x v3 and the q = states illustrated in Fig. [TJ (and 
see also, e.g., Refs. for further details). Then we 

perform a rotation of the local axes of each of the spins 
such that all spins in the reference state align along the 



2 



negative z axis. In this new set of local spin coordinates 
a complete set of mutually commuting multispin creation 
operators Cf = (CJ)^ related to this reference state is 
defined by 



|*) = 1444-..); C+ = a}, 



(2) 



where sj 



is* 



the indices n,m,k, . . . denote arbi- 



trary lattice sites, and the components of the spin oper- 
ators arc defined in the local rotated coordinate frames. 
Note that for spins of quantum number s, each site index 
in each configuration index / in Eq. ^ can be repeated 
up to a maximum of 2s times. With the set {|$),C/} 
thus defined, the CCM parametrizations of the ket and 
bra GS eigenvectors \ty) and (^| of the spin system are 
given by 



(*| = (<S>\Se~ s , 5=1 + ^ aiCj 

J#0 



(3) 
(4) 



We have defined Cq = 1, and the normalization of the 
states is clearly such that (*|*) = ($1*) = ($|$) = L 
The CCM correlation operators, S and S, contain the 
correlation coefficients, a/ and aj, which can be deter- 
mined by the CCM ket-state and bra-state equations 

i~ s He s \<l>) = ; V/ ^ 



(&\Se- s [H, C*+]e 5 |$) =0 ; V/ ^ 0. 



(5) 
(6) 



Equations ((5J and (JSJ are fully equivalent to the GS 
Schrodinger equations for the ket and bra states. They 
follow readily from the requirement that the GS energy 
functional ($|.ff|\If) be stationary with respect to varia- 
tions in all of the correlation coefficients 5/ and ai re- 
spectively (V/ 7^ 0). Each ket-state or bra-state equa- 
tion belongs to a certain configuration index I, i.e., it 
corresponds to a certain set (configuration) of lattice 
sites n,m, k, ... , as in Eq. ([3]). Using the Schrodinger 
equation, H\&) = E\ty), we can now write the GS en- 
ergy as E — ($|e _s iJe s |$). The magnetic order pa- 
rameter (sublatticc magnetization) is given by M = 
— jf Y^i=i(^\ s i where sf is expressed in the trans- 
formed coordinate system, and N(—> oo) is the number 
of lattice sites.. 

If we would be able to consider all creation and annihi- 
lation operators Cj~ and Cj , i.e. all sets (configurations) 
of lattice sites, in the CCM correlation operators S and S 
we would get, in principle, the exact eigenstate.— How- 
ever, for the many-body quantum system under consid- 
eration it is necessary to use approximation schemes in 
order to truncate the expansions of S and S in Eqs. ([3]) 
and (UJ) in a practical calculation. Then the approximate 
results for the GS energy E and the order parameter M 
will depend certainly on the choice of the reference state. 

We use for spin quantum number s = 1/2 the so-called 
LSUBn approximation scheme to truncate the expan- 
sions of S and S in Eqs. © and Q, where we include 



only n or fewer correlated spins in all configurations (or 
lattice animals in the language of graph theory) which 
span a range of no more than n contiguous lattice sites, 
where a set of sites is defined to be contiguous if every 
site has at least one other in the set as a nearest neigh- 
bor (and for more details see Refs.lHii-lii). Using an 
efficient parallelized CCM code^ we are able to solve 
the CCM equations up to LSUB10 for s = 1/2 (where, 
e.g. for the q = reference state a set of 238010 coupled 
ket-state equations has to be solved), which goes signif- 
icantly beyond earlier CCM calculations for the kagome 
HAFMjiii2& Moreover, we also use the CCM to consider 
spin quantum numbers s > 1/2. 

Since the LSUBn approximation becomes exact for 
n — > oo (as also so does the alternative SUBn-n scheme 
that we introduce and use below for values of the spin 
quantum number s > 1/2 in Sec. 1111 Bp . it is useful to 
extrapolate the 'raw' LSUBn (or SUBn-n) data to the 
limit n — > oo. There is ample experience regarding how 
one should extrapolate the GS energy per site eo(n) = 
E(n)/N and the magnetic order parameter M{n). For 
the GS energy per spin eo(n) = a + a!(l/n) 2 + a 2 (l/n) 4 
is a very well-tested extrapolation ansatzi^ ~ 38 i 40 i 46 i 48 An 
appropriate extrapolation rule for the magnetic order pa- 
rameter of highly frustrated systems i o 37 i 38 i 40 M(n) = 
bo + bjJlJnYl 2 + &2(l/rt) 3//2 - Moreover, we know from 
Refs. I37ll38ll40l that low levels of approximation conform 
poorly to these rules. Hence, we exclude the n = 2 and 
n = 3 data from the extrapolations. 

For the solution of the CCM equations we rewrite the 
Hamiltonian (fT]) in the rotated coordination frame of the 
local quantization axis 



E 

<i->i> 



f A (--fa? + 



(7) 



where we have further introduced an anisotropy param- 
eter A that now multiplies the non-Ising terms, similar 
to what was done in Refs. 04^. Note that the symbol 
< i — > j > on the sum in Eq. (J7J| now indicates a direc- 
tionality for the nearest-neighbor bonds, (and see, e.g., 
Refs. HUHI), which is different for the \/3 x -\/3 and 
the q — reference states. Starting at A = 0, where the 
corresponding reference states are eigenstates of H\, we 
can slowly increase A and hence trace the CCM solutions 
out to the true kagome point at A = 1. Moreover, A can 
be understood as a parameter that tunes the strength of 
the quantum fluctuations. 



III. RESULTS 

A. The extreme quantum case: s = 1/2 

We start with the CCM investigation of the kagome 
HAFM for spin quantum number s = 1/2. At a given 
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FIG. 2: Difference of GS energies per site, Ae = e, 
eg -0 , between the y/3 x y/3 and q = states of the spin-1/2 
kagome HAFM, for various CCM LSUBn approximations and 
spin quantum number s = 1/2. 



finite level of the CCM LSUBn scheme the treatment 
of quantum effects is performed in an approximate man- 
ner. Certainly, the treatment of quantum effects becomes 
better as the level of approximation n is increased. In 
previous studies of the GS selection based on an expan- 
sion around the classical limitiiii the -s/3 x state 
was found to be selected by quantum fluctuations. We 
present our results for the GS energy per site in Fig. [51 
where the dependence on the anisotropy parameter A of 
the difference in the energies per site between the two 



states considered, A< 



- „V3xV3 „<?=n 



, is shown for var- 



ious LSUBn approximations. Interestingly, the GS se- 
lection depends on the LSUBn truncation index n. Just 
as in linear spin-wave theory^iii there is also no GS se- 
lection (i.e., Ae = 0) at the CCM-LSUB2 level, thereby 
indicating a poor consideration of quantum effects at the 
lowest LSUBn order. As the level of approximation n 
is increased we first find that Ae < for n = 3 and 
n = 4 (i.e., the x s/3 state is selected in accordance 
with previous finding o 7 ' 11 ), but as n is further increased 
we then find that Ae > for n > 4 (i.e., the q = 
state is selected). Bearing in mind that quantum effects 
are better taken into account at higher LSUBn levels we 
might argue that strong quantum fluctuations indeed fa- 
vor the q = state. Note that this line of argument is 
also supported by our CCM results below for spin quan- 
tum numbers s > 1/2 (i.e., generally speaking, where 
quantum fluctuations are weaker), where in all levels of 
approximations the -\/3 x V3 state is selected (and see 
our discussion below in Sec. IIII B|) . 

In Fig. [3] we show the magnetic order parameter as a 
function of the anisotropy parameter A. At A = we have 
A I = s = 1/2, since the corresponding reference state is 
the exact GS of i?x=o- As A is increased the order pa- 
rameter decreases monotonically. At a certain value of 
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FIG. 3: CCM-LSUBn data for the magnetic order parameter 
M versus A for the spin-1/2 kagome HAFM, for (a) the v^x 
V 7 ^ reference state and (b) the q = reference state. For 
the extrapolations to n — > oo according to M(n) — &o + 
6i(l/n) 1/2 + b 2 {l/n) 3/2 we have used LSUBn data for n = 
4, 5, . . . , 10 as well as for n = 6, 7, . . . , 10. 



A, near the true kagome point A = 1, the extrapolated 
order parameter vanishes, thus indicating that the GS is 
magnetically disordered. The difference in the two vari- 
ants of the extrapolation (including or excluding LSUB4 
and LSUB5) may be considered as an error bar for the 
extrapolated order parameter. 

Next we analyze the model for A = 1 in more de- 
tail. The CCM-LSUBn data, as well as the extrapo- 
lated data, are listed in Table |U Moreover, we present 
results for the GS energy obtained by other methods for 
comparison. While the values eo = —0.4322 obtained 
in Ref. and eo = —0.4332 obtained in Ref. [28| can 
be considered as rigorous upper bounds of the GS en- 
ergy, the lar ge-s cale DMRG result eo = —0.4379 ob- 
tained in Ref. |28| seems to be the most accurate estimate 
presently available. The lowest extrapolated CCM en- 
ergy is eo = —0.4372 obtained for the q = reference 
state using CCM-LSUBn results for n = 6, 7, 8, 9, 10 for 
the extrapolation. This CCM estimate is very close to 
the DMRG result of Ref. HI 
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TABLE I: CCM results for the spin-1/2 HAFM on the kagome 
lattice (i.e., at A = 1). The quantity eo = E/N is the GS 
energy per spin and M is the magnetic order parameter (sub- 
lattice magnetization). The LSUBn results are extrapolated 
to n — > oo according to eo(n) = ao + ai(l/n) 2 + a2(l/n) 4 and 
M(n) = b + 6i(l/n) 1/2 + b 2 (l/n) 3/2 using LSUBn data for 
n = 4, 5, . . . , 10 as well as for n = 6, 7, . . . , 10. 



\/3 x ^3 


eo 


M 


LSUB4 


-0.408728 


0.320702 


LSUB5 


-0.414235 


0.291917 


LSUB6 


-0.418052 


0.272109 


LSUB7 


-0.420677 


0.248989 


LSUB8 


-0.423554 


0.219994 


LSUB9 


-0.424962 


0.204661 


LSUB10 


-0.426485 


0.187634 


Extrapolated (4-10) 


-0.4318 


< 


Extrapolated (6-10) 


-0.4336 


< 


9 = 


eo 


M 


LSUB4 


-0.408066 


0.322860 


LSUB5 


-0.414418 


0.286462 


LSUB6 


-0.420078 


0.248078 


LSUB7 


-0.423126 


0.225356 


LSUB8 


-0.426054 


0.202074 


LSUB9 


-0.427952 


0.186435 


LSUB10 


-0.429413 


0.172742 


Extrapolated (4-10) 


-0.4357 


< 


Extrapolated (6-10) 


-0.4372 


< 


other recent results 


Ref. 20 


-0.433 




Ref. 27 


-0.4322 




Ref. 29, N = 42 (type a) 


-0.437999 




Ref. 29, N = 42 (type b) 


-0.438143 




Ref. 28 


-0.4379 





B. Higher spin quantum numbers: s > 1/2 

Although several magnetic kagome compounds carry 
spins with s > 1/2, such as the s = 3/2 mag- 
net KCr 3 (OH) 6 (S0 4 ) 2 r£ or the s = 5/2 compound 
(H30)Fc3(OH)6(S04)2,— far fewer theoretical results are 
available for those higher spin quantum numbers. In the 
classical limit, s — > oo, thermal fluctuations may lead 
to v3 x s/3 long-range order as T — > 0&2^ In most 
papers dealing with large-spin quantum models it has 
been found that quantum fluctuations select the x 



state 



4.7.11.32 



Moreover, magnetic long-range order might 
be possible for higher spin values.—^ 

For our CCM approach for s > 1/2 we use (instead 
of the LSUBn scheme) the alternative SUBn-m approx- 
imation scheme to truncate the expansions of 5* and 
S in Eqs. © and (HJ). This is because as s increases 
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FIG. 4: Difference of the GS energies per site, Ae = e, 



V3xV3_ 





eg -0 , between the x and q — states of the kagome 
HAFM, calculated for the CCM SUB8-8 approximation and 
for spin quantum numbers s = 1/2, 1, 3/2, 2, 5/2, 3. 



the number of fundamental configuration / retained at 
a given LSUBn level also increases, since each spin at 
any site i may be raised up to 2s times by its raising 
operator sf, and hence each site index i may be re- 
peated up to 2s times in the operators Cf of Eq. ([2]). 
In the SUBn-m scheme we include no more than n spin 
flips spanning a range of no more than m contiguous 
lattice sites i 47 ' 48 In what follows we consider the case 
n = m, i.e. SUBn-n, which for s = 1/2 is identical to the 
LSUBn scheme. Since the number of coupled ket-state 
equations for a certain level of SUBn-n approximation 
increases with increasing spin quantum number s, the 
highest level of approximation we can consider is SUB8- 
8 for s = 1, 3/2, 2, 5/2, 3. The maximum number of ket- 
state equations we have to take into account is 416126 
for s = 3. For the extrapolation to n — > oo we use the 
same extrapolation formulas as for s = 1/2, and consider 
the SUBn-n data for n — 4,5, 6, 7, 8. 

First we discuss the GS selection. In Fig. [4] we present 
the dependence on the anisotropy parameter A of the en- 
ergy difference per site, Ae = e^ x ^ — e^ -0 , for the 
highest level of approximation that we have performed 
(viz., SUB8-8), and for values of the spin quantum num- 
ber s = 1/2, 1, . . . , 3. We find that the V3 x \/3 state is 
selected for all values s > 1/2. This is in agreement with 
previous studies based on an expansion around the clas- 
sical limit j^iLii but it is in contrast to our findings for the 
extreme quantum case s = 1/2. Hence, interestingly, our 
results suggest that for the frustrated quantum spin sys- 
tem under consideration the s = 1/2 case and the cases 
s > 1/2 may exhibit different behavior. It is interesting 
to note that a similar effect has also been observed for 
a frustrated quantum spin chain.— We see clearly that 
the energy difference Ae scaled by s 2 decreases mono- 
tonically with increasing s, thereby demonstrating that 
for s — > oo the V3 x y3 and the q = states become 
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FIG. 5: Extrapolated magnetic order parameter, M/s, versus 
A for (a) the x y/3 reference state and (b) the q — 
reference state of the kagome HAFM, for various values of the 
spin quantum number s. For the extrapolations to n — > oo 
according to M(n) = b + &i(l/n) 1/2 + b 2 (l/n) 3/2 we have 
used SUBn-n data for n = 4, 5, . . . , 8. [Note that even for the 
s = 1/2 case we have excluded the (available) LSUB9 and 
LSUB10 data here to be consistent with the s > 1/2 cases.] 



degenerate. 

Next we discuss the magnetic order parameter M. To 
compare results for various values of s it is useful to con- 
sider M/s. Since for s — > oo the chosen reference state 
is an eigenstate we would get M/s = 1 within our CCM 
approach in this classical limit. Hence, we may expect 
that by increasing the spin quantum number the quan- 
tity M/s becomes nonzero for a certain value s > sq in 
the whole range < A < 1 (and see also, e.g., Ref. 0). 
In Fig. [5] wc show the extrapolated magnetic order pa- 
rameter for both reference states as a function of the 
anisotropy parameter A, for values s = l/2, 5/2, 3. 
For large values of s the scaled order parameter becomes 
almost constant, M/s ~ 1, over a wide range of A val- 
ues. However, as the true kagome point at A = 1 is 
approached we find a steep decay of M/s. Nevertheless, 
only for s = 1/2 (as discussed already above) and for 
s = 1 does the extrapolated order parameter vanish at 
A = 1, whereas M/s remains nonzero for s > 1. Hence, 
our data suggest that for higher values of s a V3 x s/3 



TABLE II: CCM results for the HAFM on the kagome lattice 
(i.e., at A = 1) for spin quantum numbers 8 = 1, 3/2, 2, 5/2, 3. 
The quantity eo = E/N is the GS energy per spin and M 
is the magnetic order parameter (sublattice magnetization). 
The SUBn-n results are extrapolated to n — > oo according to 
e (n) = a +a 1 (l/n) 2 +a 2 {l/n) 4 and M(n) = 6 +&i(l/n) 1/2 + 
6 2 (l/n) 3/2 using SUBn-n data for n = 4, 5, 6, 7, 8. 





V3 x V3 


q = 


s = 1 


e /s 2 M/s 


eo/s 2 M/s 


SUB8-8 
extr4-8 


-1.383644 0.580079 
-1.4031 < 


-1.379680 0.607293 
-1.3965 < 


s = 3/2 


e /s 2 M/s 


eo/s 2 M/s 


SUB8-8 


-1.257354 0.690229 


-1.254588 0.709167 


extr4-8 


-1.2680 0.0744 


-1.2643 0.2438 


s = 2 


e /s 2 M/s 


e /s 2 M/s 


SUB8-8 


-l.iy5442 0.735642 


-1.193145 0.754580 


extr4-8 


-1.2026 0.2029 


-1.2000 0.3645 


s = 5/2 


e /s 2 M/s 


e /s 2 M/s 


SUB8-8 


-1.157697 0.766290 


-1.155703 0.785822 


extr4-8 


-1.1627 0.2942 


-1.1607 0.4586 


s = 3 


e /s 2 M/s 


eo/s 2 M/s 


SUB8-8 


-1.132263 0.788722 


-1.130497 0.808862 


extr4-8 


-1.1360 0.3583 


-1.1344 0.5256 


s — ► oo 


e /s 2 M/s 


e /s 2 M/s 


exact 


-1 1 


-1 1 



magnetic order might be possible. To provide more de- 
tailed information on the GS properties of the kagome 
HAFM we present in Table |H]CCM-SUB8-8 data as well 
as extrapolated data at A = 1 , for values of the spin quan- 
tum number s = 1,3/2,3,5/2,3. We see clearly that 
the scaled GS energy per spin approaches the classical 
value, eo/s 2 = — 1, quite rapidly as s is increased. On 
the other hand, even for the largest spin considered here 
(viz., s = 3) the extrapolated order parameter remains 
relatively small, particularly for the y/3 x y/3 state. 

Our CCM data for eo and M, available up to s = 3, 
together with the known results in the classical limit, 
liiTis-yco eo/s" = —1 and lim^oo M/s = 1, also allow 
us to discuss the s-dependence of e /s 2 and M/s in the 
large-s limit. In spin-wave theories one typically ob- 
tains expansions for e D and M in powers of l/si 53 ' 54 
For the GS energy of the kagome HAFM the standard 
linear spin-wave theory yields for both the s/3 x y/3 
and the a = states, eo/s 2 = — 1 — 0.4412/s (and see 
Refs. 00) ■ On the other hand, due to the presence 
of the flat zero mode in the kagome HAFM, the inte- 
gral for the order parameter diverges Using an ef- 
fective spin- wave theory, in which short-wavelength fluc- 
tuations are neglected, Asakawa and Suzuki^ obtained 
M/s = 1 — 0.336/s, whereas Chubukm«£ found fluctu- 
ation corrections proportional to s -2 / 3 (in contrast to 



6 



fit using s=3,5/2,2,3/2 
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FIG. 6: Dependence on the spin quantum number s of (a) 
the extrapolated scaled GS energy per spin, eo/s 2 , and (b) 
the extrapolated scaled magnetic order parameter, M/s, of 
the kagome HAFM. For M/s the exponent a — 1/2 for the 
\/3 x \/3 state and a = 2/3 for the q = state, (symbols: 
CCM data points, lines fits to the data points). 



conventional spin-waves) by using a self-consistent spin- 
wave approach. 

Here we use our extrapolated CCM data for s = 
3/2,2,5/2 and 3 to find the leading corrections to the 
classical values. By fitting the extrapolated GS energy 
eo/s 2 with a fitting function f(s) — — 1 + a\s~ a2 we 
obtain a value for the exponent a 2 very close to one 
for both reference states. Hence, finally we have fit- 
ted the extrapolated CCM data for eo/s 2 using /(,s) = 
— 1 + X1.S" 1 + x 2 s~ 2 . The fits yield the corresponding val- 
ues xi = -0.414 (-0.410) and x 2 = 0.018 (0.020) for the 
\/3 x \/3 (q = 0) reference states, as shown in Fig. EJa). 
Obviously, the | prefactor is close to that of the linear 
spin-wave theory, and the contribution of the next order 
is small. Note, however, that if the fitting function /(s), 



with the above given values for x\ and x 2 , is applied to 
the s = 1/2 case, it does not reproduce the GS selection 
of the q — state in this extreme quantum case [c.f., 
Fig-IHa) at the value 1/s = 2]. 

Next we use g(s) — l + &is~ b2 as a fitting function to fit 
the extrapolated CCM order parameters M/s, again us- 
ing data for s = 3/2, 2, 5/2 and 3. For the exponent b 2 we 
get the values b 2 = 0.529 for the -\/3 x -\/3 reference state 
and 62 = 0.666 for the q = reference state. Clearly, 



the leading correction is not proportional to 



rather 



Chubukov's result^ of a leading correction proportional 
to ,s~ 2 / 3 is confirmed for the q — state. By contrast, 
for the \/~?> x y/3 reference state our results are in favor 
of a leading correction for the order parameter propor- 
tional to s~ x / 2 . Hence, finally we have fitted the ex- 
trapolated CCM data for Mj s using the fitting functions 
g(s) = l + yis- 1 / 2 +y 2 s- 1 [g(s) = 1 + yis- 2 ^ 3 + y 2 s~ 4/3 ] 
for the \/3 x V3 [q = 0] reference states. The fits 
yield the corresponding values y\ = —1.058 [—1.000] and 
y 2 = -0.094 [0.006], as shown in Fig. E^b). 



IV. SUMMARY 

In the present investigation we present data for the GS 
energy per spin, en, and the order parameter (sublatticc 
magnetization), M, of the kagome HAFM for spin quan- 
tum numbers s = 1/2,1,3/2,2,5/2,3, using high-order 
CCM-SUBn-n calculations based on the \/3 x y/3 and 
the q ~ reference states. Our best estimate of the GS 
energy for the s = 1/2 case, viz., en = —0.4372, is clearly 
below rigorous upper bounds reported recently ) 27 ' 28 and 
it also agrees well with recent accurate DMRG— and 
results. We find that the GS selection by quantum 
fluctuations depends on the spin quantum number s. 
While for s = 1/2 the q = state is selected, for all 
values s > 1/2 the a/3 x state has lower energy. 
The order parameter M/s obtained by an appropriate 
extrapolation of the CCM-SUBn-n data vanishes for 
,s = 1/2 and s = 1, but we get small (but nonzero) values 
M/s > for all values of the spin quantum number 
s > 1. Using CCM data for s = 3/2,2,5/2 and 3 we 
determine also the leading quantum corrections to the 
classical values of the GS energy and the order parameter. 
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